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Abstract 

We present an inequality that gives a lower bound on the expecta- 
tion value of certain two-body interaction potentials in a general state 
on Fock space in terms of the corresponding expectation value for ther- 
mal equilibrium states of non-interacting systems and the difference in 
the free energy. This bound can be viewed as a rigorous version of first 
order perturbation theory for many-body systems at positive temper- 
ature. As an application, we give a proof of the first two terms in a 
high density (and high temperature) expansion of the free energy of 
jellium with Coulomb interactions, both in the fermionic and bosonic 
case. For bosons, our method works above the transition temperature 
(for the non-interacting gas) for Bosc-Einstcin condensation. 

1 Introduction 



Correlations play a crucial role in quantum-mechanical many-body systems. 
They result from interactions among the particles, and it is typically very 
difficult to obtain information about them in a mathematically rigorous fash- 
ion. Approximate theories are often arrived at by neglecting correlations, 
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for instance in Hartree-Fock theory for fermions. For the problem of esti- 
mating the validity of such approximations, it is necessary to estimate the 
magnitude of correlations present in the state of the interacting system. 

In UJ], Graf and Solovej present a correlation estimate which is applica- 
ble for the study of this problem at zero temperature, i.e., for systems in 
their ground states. The inequality presented there is motivated by earlier 
correlation estimates by Bach PQ and Bach et al. 2 . Roughly speaking, it 
estimates the difference of the interaction energy in a general state and the 
ground state of a non-interacting system in terms of the difference of their 
one-particle density matrices. Moreover, at least in the case of fermions, 
the one-particle density matrix can be easily controlled in terms of the total 
kinetic energy. For bosonic systems, the situation is more complicated, and 
the correlation estimate in is only applicable provided one can prove the 
existence of Bose-Einstein condensation — in general a very difficult task 
for interacting systems. 

With the aid of the correlation estimate just mentioned, Graf and Solovej 
were able to derive the first two terms in a high density expansion of the 
ground state energy of fermionic jellium (HI Thm. 2] with Coulomb inter- 
actions. High density corresponds to small coupling, and hence the result 
can be viewed as rigorous estimate of the validity of first-order perturbation 
theory for this system. 

In this paper, we present a method that is applicable to the aforemen- 
tioned problem for systems at positive temperature. Unlike the situation 
for the ground state, the knowledge of the one-particle density matrix alone 
does not yield much information about correlations present in the state. As 
an additional input one needs to know that the entropy of the state is close 
to the maximal value possible for given one-particle density matrix; this 
maximum is attained by the corresponding quasi-free state. More precisely, 
we will estimate the difference of the interaction energy of a general state 
and the thermal equilibrium state of a non-interacting system in terms of 
the relative entropy of these two states. This relative entropy is related to 
the difference in free energy. Our result applies to fermions at any temper- 
ature, and to bosons above the critical temperature (for the non-interacting 
gas) for Bose-Einstein condensation. 

Our main correlation estimate is stated in Theorem 01 in Section 03 Be- 
fore describing it in detail, we present an application of the inequality to 
(fermionic or bosonic) jellium with Coulomb interactions at positive tem- 
perature. We will derive the first two terms in a high density (and high 
temperature) expansion of the free energy. In the fermionic case, this result 
can be viewed as the positive temperature analogue of Theorem 2 in j^j. 
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Our estimate is general enough to be applicable to a wide range of pos- 
sible interparticle interactions. The two-body potential is required to be 
positive definite and, in particular, to be decomposable into characteristic 
functions of balls. In the case of the Coulomb potential, such a decompo- 
sition was first used in jSJ. The study in |7j provides a criterion for the 
possibility of such a decomposition for general radial functions, and thus 
provides many examples of interaction potentials which our method applies 
to. 

Acknowledgments. It is a pleasure to thank Elliott Lieb and Jan Philip 
Solovej for stimulating and fruitful discussions. 



2 Jellium 

Jellium is a model of a charged gas of either fermions or bosons, moving 
in a uniformly charged background. We assume that the whole system is 
neutral (in a sense to be made precise below) and contained in a (three- 
dimensional) cubic box of side length L, which we denote by A. We work in 
the grand-canonical ensemble, i.e., in the (anti-)symmetric Fock space over 
the one-particle space TL = L 2 (A;C n ). Here, n > 1 denotes the number of 
internal degrees of freedom, corresponding to particles of spin (n — l)/2. 

We denote by A the Laplacian on A with Dirichlet boundary conditions. 
We choose units such that H = 1 and 2m = 1, with m denoting the particle 
mass. For g > the background density and a > the square of the particle 
charge, the Hamiltonian on Fock space is 

H = H + aW, (2.1) 

where, in each iV-particle sector, 

N 



H = -Y,^i (2-2) 



i=l 



and 



N 

W = -J2g dy- r + y"i r + \(? \ dy x dy 2 - 

f^f J a \xi-y\ f^.\xi-Xj\ J AxA \yi-yz\ 



l 



(2-3) 

The last constant corresponds to the electrostatic energy of the background 
charge and is added to ensure the existence of a proper thermodynamic 
limit. 
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The quantity of interest is the free energy per unit volume at temperature 
T = given by 

/ F ' B (/3, g,a) = - lim J- In Tr exp[-/3iJ] . (2.4) 

Here, Tr denotes the trace either over the fermionic (F) or bosonic (B) Fock 
space. Existence of the thermodynamic limit in (|2.4j) was shown by Lieb 
and Narnhofer in ^01- There it was also shown that one would obtain the 
same result in the canonical ensemble with charge neutrality, i.e., fixing N 
to be g\A\. In particular, in our grand-canonical setting it is not necessary 
to enforce the charge neutrality N = g\A\ explicitly, it will be automatically 
satisfied (for the average particle number). 

There are three length scales in this problem; the mean particle distance 
f? -1 / 3 , the thermal wavelength /3 1 / 2 , and the inverse coupling constant a~ 1 . 
Hence, by simple scaling, 

f> B ((3,g,a) = ^V F 'W /8 ,l,«<f 1/8 ). (2-5) 

We are interested in the high density (and high temperature) asymptotics; 
more precisely, in large g for fixed (3g 2 > 3 (and fixed a). By the scaling 
property (|2.5|) . this corresponds to a limit of small coupling. 

For the statement of our main results, we will distinguish between the 
fermionic and bosonic cases. 



2.1 Fermions 

Let /(f(/3, q) denote the free energy (per unit volume) of a non-interacting 
gas of spin (n — l)/2 fermions, at inverse temperature /3 and average density 
g. It is given by 

f!(P, e ) = sup { W - -r^ / dp In (l + e -^ 2 -M)) ) . (2.6) 

The supremum in (|2,6|) is attained uniquely at some fi = $(f3,o). We 
denote the fugacity by z = for this value of fi. Note that z depends only 
on 13 q 2 ^. Let 

7 ° F(P) = z-ijifi + 1 ' (2 - ?) 

and let 7q (x) = (2tt)~ 3 J dp^Q{p)e ipx denote its inverse Fourier transform. 
Note that rvy^ (0) = g. 
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THEOREM 1 (High Density Asymptotics for Fermions). As g — > oo 

and (3 — * 0, 

f(P, g, a) = f!((3, g)-^l dx - o(g^) , (2.8) 

2 Jr3 \x\ 

with < o(£» 4 / 3 ) < C(Pg 2 ' 3 )ag^ 3 (ag~ 1 ' 3 ) 1 ^ 4:8 . Moreover, the function 
C(f3g 2 / 3 ) is uniformly bounded on compact intervals in (0, oo). 

Note that, for fixed f3g 2 / 3 (and fixed a), the first term on the right side 
of (|2.8|) is 0(g 5 / 3 ), whereas the second term is 0(g 4 ^ 3 ). Theorem ^ is the 
positive temperature analogue of Theorem 2 in 

We remark that (j2.8j) actually holds uniformly in (3g 2 l 3 for bounded 
l/(Pg 2 / 3 ), with possibly a worse exponent in the error term than the one 
given in Theorem ^ he., it is uniform as the ground state is approached. 
This can be proved by supplementing our lower bound with a bound ob- 
tained with the method in at very low temperatures. We do not give the 
details here, but refer the reader to ^H] where a similar argument was given 
in the case of a dilute Fermi gas with short-range interactions. 



2.2 Bosons 

For bosons we have to restrict our attention to temperatures bigger than 
the critical temperature (for the non-interacting gas) or, equivalently, to 
g < q c (J3) = n(4Tr/3y 3/2 Y,e>i£~ 3/2 ■ Let fo(/3,g) denote the free energy 
(per unit volume) of a non-interacting gas of spin (n — l)/2 bosons, given 

by 

ff(fi, ,) = sup {„ + ^ * h, (l - e-*'-»>) } . (2.9, 

For g < g c {(3), the supremum in (|2.9|) is attained at \i = ${j3,g) < 0. 
Denote the fugacity by z = er^ < 1 for this value of {l. Again, z depends 
only on the dimensionless quantity [3g 2 / 3 . Analogously to (|2.7|) . let 

7ob) = -zijf* _ i , (2-10) 

and let 7 B (x) = (2vr)- 3 / dpj}f(p)e ipx denote its inverse Fourier transform. 

THEOREM 2 (High Density Asymptotics for Bosons). .4s g — > oo 

and (3 -» (with 0g 2 / 3 < pg c (P) 2/3 ), 

/ B (/3, g, a) = / B (/3, g) + ^f dx - (, 4 / 3 ) , (2.11) 

2 Jr3 \x\ 
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with < o(g^ 3 ) < C(f3g 2 l 3 )ag i l z (ag V 3 ) 1 / 48 . Moreover, the function 
C((3g 2 ^ 3 ) is uniformly bounded on compact intervals in (0, Pq c (P) 2 ^) ■ 

As in the fermionic case, the first term on the right side of Q2.11JI is 
0(p 5//3 ), whereas the second term is 0(^ 4//3 ). Note that the second term 
diverges as g — > g c ((3)- This shows that (|2.11f) can not hold uniformly as g 
approaches the critical density, since f B ((3,g,a) < f B (oo,g, a) < for any 
f3 and g. At zero temperature, the leading term in the energy density as 
g — ► oo is actually 0(g 5 ^) [UEj- in particular, first order perturbation 
theory (in the grand canonical ensemble) is not applicable below the critical 
temperature, due to the large fluctuations in particle number. These large 
fluctuations cannot be present in the interacting system, for any non-zero 
value of the coupling parameter a. 

The key ingredient in the proof of Theorems Hand El is a new correlation 
estimate, which we present next. 

3 Correlation Estimate 

In this section, we will describe our main correlation estimate, which will 
then be used in the proof of Theorems ^ an d EJ For £ £ M 3 and r > 0, 
let Xr,§ denote the characteristic function of a ball of radius r centered 
at £. The function Xr£ defines a projection operator on L 2 (IR 3 ;C n ) and 
also, in a natural way, on the subspace L 2 (A;C n ). Let n r ^ denote the 
operator on Fock space that counts the number of particles in this ball, 
i.e., the second quantization of the projection Xr£ on TL = L 2 (A;C n ). Our 
correlation estimate concerns a lower bound on the expectation value of the 
number of pairs of particles inside a ball of radius r or, more precisely, on 

/ d£TrK,cK $ -±)r] . (3.1) 

Here, T is a density matrix, i.e., a positive operator on Fock space with trace 
equal to one, defining the state of the system. 

Let 70 denote the one-particle density matrix of a (grand-canonical) non- 
interacting (Fermi or Bose) gas at temperature T = with chemical 
potential \i = /Xq' B (/3, g), as defined after Eqs. 1)2. 6j) and (|2.9jl . respectively. 
We choose periodic boundary conditions for 70, which has the advantage 
of 70 having a constant density. Note that the choice of \x implies that 
|A|f? = tr7o = |A|g + o(|A|) in the thermodynamic limit. Here and in the 
following, we denote the trace over the one-particle space Ji = L 2 {K\ C n ) by 
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tr, whereas the trace over Fock space is denoted by TV. The kernel of 70 is 
given by 

l0 ( x ,a-y,T) = ±- £ 7o F ' B (p)e^)^ T , (3.2) 
1 1 pe^z 3 

where 7q' B (p) is given in (|2.7|) and (|2.1()j) . respectively, and a and r label 
the spin states. 

Let To denote the quasi-free state on Fock space with one-particle den- 
sity matrix 70. It is the Gibbs state (at inverse temperature 1 and chemi- 
cal potential 0) for a non-interacting system with one-particle Hamiltonian 
ln[(l =p 7o)/7o]- Here and in the following, =p means — for fermions and + 
for bosons (and vice versa for ±). We note that for r = Fq, the expression 
in ()3.1|) can be easily calculated. Namely, for any r and £, 

Tr K^(n n5 - l)r ] = (tr [Xr,£7o]) 2 T tr (Xr,§7o) 2 • (3.3) 
Hence, after integration over £, 

/ d£ Tr [ra r £ (n r ^ - 1) r ] 

= dxdyJ r (x-y)[Q 2 TT,a\'yo(x,cr;y,a)\ 2 ] , (3.4) 
JAxA 

where we denoted J r (x) = J dy Xr,^(y)Xr,^(x — y)- (Note that J r is inde- 
pendent of £.) Here, we have also used that 7 o has a constant density 
Q = Ea 7o(af, cr; o") for x G A. 

We want to show that for states T that are in some sense close to the 
state To, the expectation value (|3.1j) is close to (|3.4|) . A convenient way to 
characterize this "proximity" is the relative entropy: For two general states 
r and T on Fock space, the relative entropy is given by 

S(r,T) = Trr(lnr-lnT). (3.5) 

Note that < S(T, T) < 00. Although S does not define a metric, it 
measures the difference between two states in a certain sense. In particular, 
S dominates the trace norm. More precisely, S(T,T) > 2||T — T||^ |14l 
Thm. 1.15]. 

Note that the relative entropy can also be interpreted as a difference in 
free energies. More precisely, if T = exp(— (5{H — F)) for some (3 > 0, with 
F = — f3~ l In Tr exp(— (3H) the corresponding "free energy", then 

/9 -1 5(r, T) = Tr [HT] + ^ _1 Tr r In T - F . (3.6) 
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Note that — TrTlnr is just the von-Neumann entropy of T. Hence the 
first two terms on the right side of (|3.6j) correspond to the free energy of 
r (with Hamiltonian and temperature determined by T), whereas F is the 
free energy of T. 

Our main result estimates the difference of the expectation value 1|3.1|) 
for Tq and a general state T in terms of the relative entropy S(T, Tq). More 
precisely, the following Theorem, which is the main new result of this work, 
holds. 

THEOREM 3 (Main Correlation Estimate). Let To be given as above, 
with one-particle density matrix 70 and density g, and with fj, £ R for 
fermions and \i < for bosons. Let T be any other state on (fermionic 
or bosonic) Fock space. For any 2r < d < L/2, we have that 



Here, C z ' are constants depending only on z = e"^, which are uniformly 
bounded on compact intervals in (0, 00) and (0,1), respectively. 

We emphasize again that, according to (|3.4[) . the second line in (|3.7|) 
equals the first in the case T = Tq. Although the inequality (|3.7j) is not 
sharp in this case, the parameter d can be made very large to obtain an 
error with is, in the thermodynamic limit, of lower order than the volume. 
(The restriction d < L/2 in Theorem 01 is purely technical and could in 
principle be avoided by a slight modification of the proof. Since we are 
mainly concerned here with the application of (j3.7j) in the thermodynamic 
limit L —* 00, we have refrained from doing so.) 

Note that Theorem |31 gives an estimate on a "local" quantity, like the 
expectation value of the number of pairs of particles inside a small ball, in 
terms of a "global" quantity as the relative entropy. The strong subadditivity 
of entropy plays a crucial role in this estimate. Before we give the proof of 
Theorem|21 we show how it can be used to prove the applications to Coulomb 
systems stated in Theorems ^ and El 




(3.7) 



(1 + r 3 g) | A| 3 / 4 d 3 (l + (3d~ 2 ) S(T, T ) + ^d^A] 



1/4 



S 



4 Proof of Theorems [T] and [2] 



We are going to treat the fermionic and bosonic case simultaneously, merely 
pointing out the differences if necessary. We start by deriving a lower bound 
on the free energy. Note that if T denotes the Gibbs state of H at temper- 
ature (and zero chemical potential), then charge neutrality (as proved 
in \1U\ ) implies that 

lim — TrAT = p (4.1) 

L-*oo |A| 

for any fixed (3 and a > 0. Here, N denotes the number operator on Fock 
space. Application of the Peierls-Bogoliubov inequality then leads to the 
lower bound 

f' B (P, g,a) > f^ B ((3, g) + alimsnp-^TrWT . (4.2) 

L^oo | A | 

To estimate the expectation value of W in the Gibbs state T, we will split 
the Coulomb potential into a long and short-range part. 



4.1 Long- Range Part 

We write the Coulomb potential as [3] 

1 1 f°° 1 f 

f - v\ n Jo r Jr 3 

As in Section Xr,£ denotes the characteristic function of a ball of radius r 
centered at £ € M 3 . We split the r-integration into a part r < R and a part 
r > R and, correspondingly, write 

1 



F - y\ 



V <R {x -y) + V >R (x - y) . (4.4) 



Note that V < r{x) = for |x| > 2R. For the long-range part V>r, we note 
that it has a positive Fourier transform, as follows immediately from the 
decomposition (|4.3|) . Hence we obtain the lower bound ^3 4.5.20] 



^2 V >R (xi -Xj)>^2g / dy V >R (xi - y) (4.5) 

' A ^ AT A — 1 ''A 

-,g 2 I dyi dy 2 V >R {y 1 - y 2 ) - ttV >r (0) . 
JAxA 1 



l<i<j<N i=l 



This estimate actually holds for any g > 0. The last term equals V >R (0) 
4/(3R), and hence will be negligible if we choose R 3> g^ 1 ^ . 
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4.2 Short-Range Part 

As in Section let n r ^ denote the operator that counts the number of 
particles in a ball of radius r centered at £, i.e., the second quantization of 
the projection Xr£ on Tt = L 2 (IR 3 ; C n ). The expectation value of the short- 
range part V < r of the interparticle interaction in a state T on Fock space 
can be written as 

-!- f R dr\ [ d£ Tr K, e (n r , e - 1) T] . (4.6) 
^ JO r JR 3 

For a lower bound, we can now apply our main correlation estimate, Theo- 
rem EH to the expression (|4.6|) for any fixed r. 

Recall that 70 denotes the one-particle density matrix of a non-inter- 
acting (Fermi or Bose) gas at inverse temperature 0, with chemical poten- 
tial fi = /Xq' B (/?, g), and with periodic boundary conditions; To denotes the 
corresponding quasi-free state on Fock space. Theorem El states that for any 
2r <d< L/2, 

/ d£ Tr [n r ,£ (n r £ - 1) V] 
Jr 3 

> / dxdyJ r (x - y) [g 2 =F 
JAxA 

-Cj>V£(l + r 3 ^)|A| 3 / 4 

For r the Gibbs state of H, an upper bound on S(T, Tq) is, in fact, easy 
to obtain. Using the fact that the quadratic form domain of the Dirichlet 
Laplacian is contained in the quadratic form domain of the Laplacian on A 
with periodic boundary conditions, we can write 

5(r,r ) =/?Tr (H - fiN)T + Tr T In T =F tr ln(l T To) 

= - In Tr exp[-0H] - /3fx Tr NT - Pa Tr WT =F tr ln(l =F To) • 

(4.8) 

We now use the lower bound W > — const. A| f? 1 / 3 [Hlj, as well as the 
fact that |A| -1 TriVT — > g in the thermodynamic limit, as explained in the 
beginning of this section. This leads to the estimate 

S(T,T ) < \A\P^ B (P,g,a)-f^ B (P,g)) + const. /?|A|a// 3 + o(|A|) . 

(4.9) 



Y,aho( x i a 'iyi a )\ 2 ] 

> (x + 0d -2} S (T, T ) + /3V2 d -i| A | 



(4.7) 

1/4 
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As the upper bound to the free energy in Section [4.41 shows, the first term 
is negative in the fermionic case and can thus be neglected for an upper 
bound. In the bosonic case, it is bounded above by C z f3\A\ag^ 3 for some 
constant depending only on z. This follows immediately from the upper 
bound leading to (|2.11|) . together with simple scaling. (Note that C z diverges 
as z — > 1). Hence, in general, 

5(r,r ) < C z f3\A\ag A / 3 + o(\A\) . (4.10) 

Here and in the following, we abuse the notation slightly and denote by 
C z any expression that depends only on z (and is uniformly bounded on 
compact intervals in (0, oo) in the fermionic case and (0, 1) in the bosonic 
case). 

We insert the bound P~l"0|) into (|4~7|) . Choosing d = l3~ 1/% a~ 1/A g~ 1/3 
we thus obtain that, as long as d > 2r (and ag^ 1 ^ 3 < const. ), 

I d^^Kf Kc-i)r] 

J®? 

> / dxdyJ r (x-y)[g 2 ^^2Jj (x,a;y,a)\ 2 ] 

J AxA 

-C z r 3 g(l + r 3 g) \A\ (/? 5 / 2 a/ /3 ) ^ + o(| A\) . (4.11) 

Here, we have also used that g = g + o(l) as L —* oo. Note that we are 
going to use this estimate in (|4.6|) only for r < R. Below we will choose 
R <. /3 _1 / 8 a -1 / 4 ^ 1 / 3 , hence (j4*?TT|l will be applicable. 

For a lower bound, we can restrict the r-integration in (|4,B|) to Rq < 
r < R for some < Ro < R, and simply neglect the contribution from 
the r < Rq part. A simple estimate, using that J2 (T \'yo{x,a;y,a)\ 2 < g 2 , 
shows that in the state Tq this r < Rq contribution is bounded above by 
(27r)- 1 |A|(47r/3) 2 ^ 2 J R§. We then have 

-!- / dr -1 / di Tr [n rA (n r ^ - 1) T] 
27 r Jo r Jm* 

> \ f dxdyV <R (x-y)[g 2 TEMx,a ]y ,a)\ 2 ]-\A\^g 2 R 2 
1 J AxA y 

- C z g\A\ (±- + R 2 g^j [ag- 1 / 3 ) 1 ^ + o(|A|) . (4.12) 
Here, we have used again that z is a function of f3g 2 ^ 3 . 
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4.3 Final Lower Bound 

For the one-body part containing V<r (i.e., the interaction with the back- 
ground), we can use the simple lower bound 

- V e I dx V <R (x - Xi ) > -Nq I dx V <R {x) . (4.13) 

~l Ja Jr 3 

In combination, ()4.5|) . (|4.12|) and (|4.13|) yield, in the thermodynamic limit, 

hminf-Lwr> T - I dxV <R {^)\ll ,B ^)\ 2 - — -—Q 2 Rl 
l^oo |A| " 2 j M 3 < K JU0 Wl 3R 9^ 

-C zg (±- + R* g y ae -V 3 ) 1/m . (4.14) 

In the fermionic case, we can simply use V <R {x) < for a lower bound. 
In the bosonic case, we write V <R (x) = |a?| 1 — V >R (x) and estimate (using 
V >R (x) < V >R (0) = 4/(3i?)) 

5/ dxy>K(x)|7 B (x)| 2 <||(2vrr 3 / dp | 7o B (p)| 2 < ||-L- . (4.15) 
2 Jr3 ait Jr3 6KI - z 

With the choice # = g- 1 / 3 (ag~ 1 / 3 y 1 / AS and i? = l/(R 2 g) this yields 



liminf ^TVWT > T- / ' ,u , y ' - C 2 ^ d (4.16) 
i-»oo |A| 2 J r3 \x\ 

for some constant C z depending only on z. Inserting this bound into (|4.2|) 
finishes the proof of the lower bound. 

4.4 Upper Bound 

For the upper bound to the free energy, we use the variational principle, 
which states that 

-4hiTr exp[-pB\ < Tr HT - ^S(T) (4.17) 

for any state T on Fock space. Here, S(T) = — TrTlnT denotes the von- 
Neumann entropy. We choose as a trial state T a quasi-free state with 
one-particle density matrix 7 given by the kernel 

v, V, t) = g{x)g{yYyl' B (x - y)8 (TjT . (4.18) 
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Here, < g(x) < 1 is a continuously differentiable function with the property 
that g(x) = for x A, g(x) = 1 if x € A and dist(x,<9A) > R, and 
|V<7| < const. R . We shall choose the variational parameter R to satisfy 
l/(Lg 2 / 3 ) < R < (L^ 2 )- 1 / 5 for large L. 

The calculation of the energy of the state T is similar to the correspond- 
ing calculation in 6 . It is in fact simpler since the particle number does not 
have to be fixed. 

Let ip(x) = 7(x, a; x, a) denote the density of 7. A simple computa- 
tion (compare with (|3.3[) - (j3,4[) ). using the fact that T is a quasi-free state, 
yields 

Triyr = I f dxdy-^— [(<p(x) - e )(<p{y) - q)t V^)\ 2 } • 

jaxa \ x y\ 

(4.19) 

Note that (p(x) = g(l—g(x) 2 ) and hence, by definition, ip(x) = g if x £ A and 
x is at least a distance R away from the boundary of A. Using the Hardy- 
Littlewood-Sobolev inequality _9, Thm. 4.3], it is easy to see that the first 
term on the right side of (|4,19|) is bounded from above by const. g 2 (L 2 R) 5 ^ 3 
and is thus negligible in the thermodynamic limit, if R <C (Lg 2 )^ 1 /^ . In the 
fermionic case, the second term is bounded from above by 

~ I dxd yT ^—\^{x-y)\ 2 (l-2{l-g{x) 2 )) , (4.20) 

* jaxa \ x y\ 

which yields the desired expression in the thermodynamic limit, provided 
R <C L, which is amply satisfied for our choice of R. In the bosonic case, 
we can simply use g < 1 to obtain the desired bound. 
The kinetic energy of T is given by 

Tr# r = -nA7o' B (0) / dxg(x) 2 + g[ dx\Vg(x)\ 2 

nT 2 

< -n|A|A7^' B (0) + const. ^— . (4.21) 

R 

Again, the first term is the desired expression, and the last term is negligible 
if i?> l/(Lg 2 / 3 ). 

It remains to derive a lower bound on the entropy S(T). We claim that 

x J 7o ,B (p) ln7o ' B (P) ± (l T 1q' B (p)) In (l T 7o ,B (p) 

(4.22) 
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which gives the desired quantity as long as R <C L. Inequality 1)4.22(1 follows 
from a variant of the Berezin-Lieb inequality [3J El ■ The one-particle density 
matrix (|4.18f) can be written as 

7 = y2 dpio' B (p)g\p,v)(p,v\g , (4.23) 

where g denotes multiplication by g(x), and \p, a) denotes a plane wave with 
wave function (27r) -3 / 2 exp(ipx) and spin a. Moreover, since T is a quasi-free 
state, 

5(r)=tr S ( 7 ), (4.24) 

where we denoted s(t) = — tint =p (1 =p t)ln(l =F t) for t > 0. Note that s 
is a concave function, with s(0) = 0. Hence we can apply the Berezin-Lieb 
inequality, in the form proved in Thm. Al in 16 . Noting that 



Y, f dpg\p,a)(p,a\g = g 2 <1, (4.25) 



as well as (p, cr\g 2 \p, a) = (2tt) 3 J dxg(x) 2 , this yields (|4.22|) . 
We conclude that 



- lim _L lnTr exp[-/3fll < =p— / dx l7 °', ^ - nA7 n F ' B (0) 
l^oo (3\A\ 1 2 7 R3 |x| 

+ ?2^3g j£ s ^ ^' B (P) ln 7o ' B W ± (l T 7o ' B (^)) ln (l T 1o' B (p) 

(4.26) 

The last two terms together are just /J' B (/9, g). We have thus established 
the desired upper bounds. This concludes the proof of Theorems ^ and [21 



5 Proof of Theorem [3] 

5.1 Localization of Relative Entropy 

If X denotes a projection on the one-particle space H., then states on the 
Fock space can be restricted to the Fock space over the subspace XTt of 
7i. We denote such a restriction of a state F by Fx- Since Xr,£ defines a 
projection on 7i = L 2 (A; C"), we can write 

Tr [n ng (n r ,£ - 1) T] = Tr [n n€ (n n€ - 1) F Xr ( ) , (5.1) 
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the latter trace being over the Fock space over Xr,(7~L- 

It is well known ^1] that the relative entropy decreases under restriction. 
More precisely, for any two states T and T on Fock space, 



S(T,T)>S(T X ,T X ). 



(5.2) 



This property is closely related to the strong subadditivity of the von- 
Neumann entropy 113j . 

Let rj : R 3 i-> R be a function with the following properties: 

• 77 € C 4 (M 3 ) 

• r/(0) = 1, and rj(x) = for |x| > 1 

• v(p) = / dxr](x)e~ tpx > for all p £ R 3 . 

We note that such a function (with any degree of regularity) can, for in- 
stance, be obtained by taking a smooth function of compact support, and 
convolving it with itself. The resulting function is then smooth, has com- 
pact support and positive Fourier transform. In our application, we need 
the existence of the fourth derivatives at the origin (see Eq. (|5.3U|) below) . 
Given such a function 77, we define %(x) = rj(x/d) and 



Note that rfa is a periodic function with period L and, since L > 2d by 
assumption, we have that T]^ er < 1. Moreover, we define a one-particle 
density matrix 7^ on H by the kernel 



with 70 defined in l|3.2|) . This defines a positive operator, with plane waves 
as eigenf unctions, and eigenvalues determined by the convolution of % and 



If [L/2d] denotes the largest integer < L/2d, define J by L/2d = [L/2d]. 
Then d < d < 2d. For < r < d/2, let X r denote the characteristic function 
of a collection of balls of radius r, separated by 2d: 




(5.3) 



ld{x, cf\ y, t) = ^{x,a;y,T)rf d ev {x - y) 



(5.4) 







(5.5) 



^g2<iZ 3 n[0,L) 3 



where we denoted 




(5.6) 
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Note that the minimal distance between the balls is 2d — 2r > d. Hence 

X rld X r = £ x£7ix£, (5.7) 

^G2dZ3n[0,L) 3 

the off-diagonal terms vanish since 17,2(2;) = for |x| > d. I.e., X r ^ d X r is a 
direct sum of one-particle density matrices on x^^TL for £ G 2cffi 3 n [0, L) 3 . 

Let denote the quasi-free state on Fock space with one-particle density 
matrix 7^, and let V denote any other state on Fock space. The characteristic 
function X r defines a projection operator on the one-particle space Ti = 
L 2 (A;C n ). Hence the monotonicity of the relative entropy implies 

S(T,r d )>S(T Xr ,T d>Xr ), (5.8) 

where Tx r and T d; x T denote the states restricted to the Fock space over 
X r TC, respectively. Note that the one-particle density matrix of the quasi- 
free state T d: x r is given by X r ^ d X r . Hence (j5.7j) shows that T dt x r can be 
written as a product of states on the Fock spaces over the one-particle spaces 
Xr^'H- f° r £ ^ 2dl? n [0,L) 3 . Under this condition S is superadditive, as 
follows easily from subadditivity of the von-Neumann entropy [Tl|- More 
precisely, 

s(v,r d ) > s(T Xr ,r d! x r ) > £ s(r x ^,r d<i ) . (5.9) 

^G2dZ3n[0,L) 3 

We can repeat the argument above with a projector defined by the mul- 
tiplication operator X r (x + a) for some vector a S [0, 2cZ] 3 . Averaging over 
a then yields 



S(T,T d )>—j- da V Per ,T d per ) 

= w/ A ^ 5(r ^' r ^ } - (5 - 10) 



Remark. We emphasize that in order to obtain the superadditivity of 
the relative entropy leading to (|5.9|) . we have used the fact that T dj x r = 

r djX per. Our estimate applies to any density matrix having this prop- 
erty. For a general state, however, it will be difficult to check this property; 
in the case of a quasi-free state considered here, it simply translates to the 
vanishing of off-diagonal terms in the one-particle density matrix (more pre- 
cisely, the validity of (|5.7jl ). 
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5.2 Upper Bound on Relative Entropy with Cutoff 



In the previous subsection, we have shown how to localize relative entropy in 
the case when the second argument is a state that has been cut off in such 
away as to avoid correlations between balls of a certain distance. In the 
following, we will quantify the effect of this cut-off on the relative entropy. 

If T 7 denotes the quasi-free state with one-particle density matrix 7, 
and T is any other state on Fock space, then S"(T,r 7 ) is convex in 7. This 
follows from operator-concavity of the logarithm and S(T, T 7 ) = Tr Y In T — 
tiui hi7=Ftr (l=Fw) ln(l=P7), where u> denotes the one-particle density matrix 
of T. Note that 7^ can be written as a convex combination of the form 

7<i=TT7 Yl Vd{q)l^q 

2TT rjj 



|A| 

where 70,5 is defined by the kernel 

( v 1 ^ 1 



|A| 



E 



+ 7 IP- 



Hence convexity implies that, for any state T, 



S(T,r d )<±- ]T r} d (q)S(T,T 0jq ), 



(5.11) 

Jp( x -y)g cTT . 

(5.12) 

(5.13) 



qe 1 



where Tq^ denotes the quasi-free state corresponding to the one-particle 
density matrix jQ tq . 

Recall that Tq denotes the quasi-free state on Fock space with one- 
particle density matrix 70 given in (J3.2|) . i.e., Tq = Fq^. We claim that, 
for any t > 0, 

1 1 



S(r, r 0>g ) < (1 + 1- 1 ) 5(r, r ) + tr (h q - h 



where h q = ln[(l =F 7o,g)/7o,g]- 



e (l+t)ho—th q _|_ 1 e h q _|_ 1 

(5.14) 

In the Bose case, we have to assume that 



(1 + t)ho — th q > 0, which is satisfied for t small enough, as our estimates in 
Lemma|21below will show. Inequality Q5.14|) follows from the two inequalities 
(where 7 denotes the one-particle density matrix of T) 

tT-y((l + t)h -th q ) +Trrinr (5.15) 

> Ttr In (l ± e -0-+t)ho+K 



> =Ftr In ^1 ± e- ho ) + t tr (h - h q ) [e 



(l+t)ho-thq 
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and 

=Ftr In (l ± e~ h A > =Ftr In (l ± e' ho ) + tr (h q - h )[e h " ± l]" 1 . (5.16) 



Dividing (|5.15f) by t and adding (|5,16l) yields (|5.14|1 . 

To estimate the last term in (|5.14|1 . we need the following simple lemmas, 
estimating the expression 

W=>-^i^P^. (5.17) 
7o (P + Q) + 7 (P - Q) 

Note that h F ' B (p) = P(p 2 ~ /-«)■ 

Lemma 1 (Fermions). Let D z = sup u>0 [zu j (e u + z)\. Then 

-2f3q 2 (3D z + 2(3p 2 ) < h F q (p) - h F (p) < 2(3q 2 (l + 2D Z ) . (5.18) 
Moreover, 

P(q 2 ~ 2\pq\) < h F q (p) - h F (p) < (3(q 2 + 2\pq\) (5.19) 
independently of z. 

Lemma 2 (Bosons). Let D z = swp u>0 [z 2 ue u / (e u — z) 2 ]. Then 

-2pq 2 (3D Z + 2(3p 2 ) < h B (p) - h B (p) < (3q 2 . (5.20) 

Moreover, 
independently of z. 



h B (p)-h B (p)>(3(q 2 -2\pq\) (5.21) 



We defer the proof of Lemmas ^ and |2] to the appendix. 
The last term in (|5. 14f) is given by 



n 



E (^ b (p)-^' b (p) w 1 



e (l+t)hg' B (p)-thg' B {p) _|_ i e hq' B (p) _|_ l 



(5.22) 

A simple estimate on the derivative of the last term in brackets with respect 
to h F ' B (p) — h F,B (p) shows that (|5.22j) is bounded above by 



(l + t)C z E (h F q ' B (p)-h F > B (p)Y sup ' 



l<S<t e ( 1 + S ) /l o' (P)- Sh q' (?) ± 1 

(5.23) 
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where C z = 1 for fermions and C z = (1 — z) 1 for bosons. (Here we have 
used that 1 + jfi(p) < (1 - z) -1 -) The upper bounds in (jS~T%() and ()5~2U)) 

show that, for < s < t, 

(i + s )^' B (p) - sh^( P ) > { %<r\ ~ fffi 1 + 20 J ! or J*™ > 

v ' u v y 9 v y [ hQ{p)-tl3q z for bosons. 

(5.24) 

We choose t = min{l, (2/3g 2 (l + 2D Z )) 1 } in the fermionic case, and t = 
min{l, — /i/(2g 2 )} in the bosonic case. With this choice, 1)5.24(1 becomes 

(i + .)cm - > { : ^ 1 £ :™ s • <^> 

for < s < t. For — 1 < s < we use the lower bounds in (|5.19j) and 1)5.21(1 . 
respectively. It is then easy to see that in this case 

(1 + s)h F > B (p) - sh F q > B (p) > (3 [mm{p 2 , (p - q) 2 , (p + q) 2 } - M ] . (5.26) 

Applying the bounds ()5.26() and 1)5.25(1 to the denominator in 1)5.23(1 and 
using 1)5.18)) and 1)5.20(1 . respectively, to bound the expression (h F,B (p) — 
hq (p)) 2 from above, we obtain that 

Km < C z \A\(3 1/2 q 4 (5.27) 

as long as j3q 2 < const. Here we have also used that t < 1 by definition. 
(Again, as in Section 21 we abuse the notation slightly and denote by C z 
any expression that depends only on z.) 

It remains to show that 1)5.27)1 holds also for large values of (5q 2 . To do 
this, we can go back to 1)5.22(1 and apply the bounds above directly to this 
term. In case h F ' B (p) > h F,B (p), we use 1)5.25(1 (with s = t) as well as the 
upper bounds in 1)5.18(1 and 1)5.20)1 . For the case h q ' (p) < h ' (p), we use 
l)5.2fi)l and the lower bounds in ()5.19)) and (j5.21)) . We then split the sum 
into three regions according to where the minimum in 1)5.26(1 is attained, 
and change variables from p to p — q or p + q, respectively. In this way we 
see that 

22)<C z |A|£- 1 M(l + /3 1 / 2 |gf) (5.28) 



for any value of q. Hence, in particular, l)5.27|) holds for all q. 
We have thus shown that 

5(r, r , 3 ) < 2 (i + c' z p q 2 ) s(r, r ) + c z \A\p l / 2 q 4 , (5.29) 
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with C' z = 1 + 2D Z for fermions and C' z = — 1/lnz for bosons. We insert 
this bound into (|5.13|) and sum over q. We can use 



(5.30) 



and similarly for g 4 replaced by g 2 . This leads to the result that, irrespective 
of whether we consider Fermi or Bose symmetry, 



i + pd- 2 ) s(r,r ) + |A|^ 



(5.31) 



S(T,T d ) <Q 

with C z a constant depending only on z = e^. 
5.3 Final Steps in the Proof 

If n r £ denotes the operator that counts the number of particles in a ball of 
radius r centered at £, we want a lower bound on the expression 



(5.32) 



For a lower bound, we can replace the positive operator n r ^(n r ^ — 1) by 
fK(ji r £(n r £ - 1)), where 



for some K > to be determined. Then 



i for t < K 
K for t > K 



(5.33) 



Tr [n r ^ (n r4 - 1) r Xr J > Tr [/# (n r , f (n r>5 - 1)) r Xp<€ ] 

> Tr [f K (nr,e(«r,$ - 1)) r d , Xr J 



Alir 



(5.34) 



Next we note that t — /i<r(t) = [t — K] + < t 2 /(4K), and hence 
Tr [f K (n r ^(n r ^ - 1)) r d|Xr . £ ] 



> Tr [n r ,c("r,e ~ 



Tr [n^(n r , 5 - l) 2 r d>Xr J . (5.35) 
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Note that F dtXr e 1S a quasi-free state. Hence (compare with (|3.3fl — (|3.4|> ) 
/ a?£Tr [n r ^(n r ^ - l)r dtXr( ] 

-J ^-.j^M^fl. 

JAxA 

Moreover, the last term in (|5,35|) is easy to estimate. Since T dtXr £ is quasi- 
free, it can be the explicitly expressed in terms of Xr£ldXr£- A simple 
estimate then yields, in the fermionic case, 

Tr [n^(n r , f - l) 2 r d , Xr J < (tr [xrtld]? (tr [x r ,{7«i] + 2) 2 • (5.37) 

In the bosonic case, we obtain 

Tr [n^(n r , ? - l) 2 T d>x J < 24 (tr [Xr^dlf (tr [x r ,{7d] + if • (5-38) 

Note that tr[x r ,§7d] = 47rr 3 ^/3 as long as A contains the ball of radius r 
centered at £, since 7^ has a constant density g. For any £ and r we have 
tr[;\>,£7d] < 47rr 3 ^/3. Integrating over £ thus yields 

/ ci£Tr [n^(n r , s - l) 2 r d , Xr J < const. |A| (r 3 ^) 2 (l + r 3 ^) 2 . (5.39) 

JR 3 

To estimate the last term in (|5.34|) . we first note that ||T XrS — r diXr( || 2 < 
25 , (r Xr . e ,rrf jXr5 ) O Thm. 1.15]. Using Schwarz's inequality for the £- 
integration yields 

/ dZ\\T Xrt( -r d>x j\ 1 <V2(L + 2r) 3 / 2 ( [ dtSF^r^)) . 

(5.40) 

Here we have also used the fact that the integrand is zero if the distance 
between £ and A is bigger than r, since there are no particles outside A and 
hence both restricted states are the Fock space vacuum in this case. To esti- 
mate the last term in (j5.40|) . we would like to use (|5.10j) . We note that, again 
by monotonicity of the relative entropy, S(T Xri ,T diXrS ) < S(T x p^ ,T d ^p^). 

The latter quantity is periodic in £, with period L. Moreover, since r < L/2 
by assumption, the cube of side length L + 2r is contained within 3 3 copies 
of A, and hence 

^ 3 # S(T XriV r diXr £ ) < 3 3 J d£ S(r x pe r , T d<l ) . (5.41) 
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Using ()5.1(J|) this yields 

/ d£||r T ||a <4(L + 2r) 3 / 2 (3d) 3 / 2 S(r,r d ) 1 / 2 . (5.42) 

Note that (L + 2r) < (3/2)|A| 1 / 3 since 2r < L/2 by assumption, as well as 
d < 2d. 

Collecting all the terms and optimizing over K, we obtain the lower 
bound 

/ ^TrK, f K, f -i)r] 

> / dxdy J r (x - y)[f ^Y,a\ld{x,a;y,(T)\ 2 ] 

iAxA 

- const. r 3 g (1 + r 3 g) | A| 3 / 4 d 3 / 4 S(I\ I^) 1 / 4 . (5.43) 

Note that |7d(x, a; y, a)\ < \ jo(x,a;y,a)\ because of f|5.4|) and the fact that 
|f?^ er | < 1- Hence (|5.43p . together with (|5.31jl . proves the Theorem in the 
fermionic case. 

In the bosonic case, we have to estimate, in addition, the term 

dxdyJ r (x-y)h (x,a-y,a)\ 2 (l-riY T (x-y) 2 ). (5.44) 

J AxA 

We use that J r (x) < (47r/3)r 3 and \^yo(x, a; y,a)\ < g/n. Moreover, we can 
estimate r]^ er (x) 2 > 1 — const. (x/d) u for any < v < 2. Choosing v = 1/4 
we obtain the bound 

(|5.44j) < const. r 3 od~ 1//4 / dx dy \jq(x, a; y, a)\ \x - y\ 1 ^ . (5.45) 

iAxA 

By simple scaling, the integral is bounded above by C^A)/? 1 / 8 for some z- 
dependent constant. Hence the error term (|5.44() can be absorbed into the 
error terms already present in Q5.43j) . merely adjusting the constant. 
This finishes the proof of Theorem El 

A Appendix 

Proof of Lemmas Q and We first prove (|5.19|) and (|5.21|) . Since both x 
ln[(2 — x)/x] and x i— > ln[(2 + x)/x] are monotone decreasing (for < x < 2 
and x > 0, respectively), we can obtain upper and lower bounds on h^' B {p) 
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by replacing 7q ,B (p + q) and 7o' B (p — q) by the minimal and maximal value 
of these two expressions, respectively. This yields (|5.19(l and 1)5.21(1 . 

The upper bound in Q5.2UJI follows immediately from convexity of the 
map x i — > ln[(2 + x) jx\ for x > 0. 

The proof of (|5.18|) and the lower bound in (|5.2U|) is a bit more tedious, 
but elementary. For convenience we set = 1, the correct /3-dependence 
follows easily by scaling. For < A < 1, we define /(A) = h x ' (p). Note 
that /'(0) = and hence 



K' a (p) - K' D ( P ) = /(l) - /(0) = f 1 d\ (1 - A)/"(A) . 

J o 

To calculate /"(A) it is useful to note that 

<?V7o' B (p) = -2p<77o' B (p)(l =F7o' B (p)) 

and 



+ 8(pq)Vo' B (p)(l T 7o F ' B W) (| -F 7o F ' B W) • 



Denoting p± = p± \q and 7 ± = 7o' B (p±), we therefore have 
1 1 



/"(A) 



+ 



(2 T7+T7 _) 2 (7+ + 7-) 



1 



1 



+ 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



x (2p+<?7+( 1= F7+) -2p_g 7 _(l=F7_ 

-2g 2 7 + (lT7+)-2 ( zV(lT7-) 



2 =F 7+ T 7- 7+ + 7- 
+ 8(p+g) 2 7+ (l =F 7+)(l =F 7+) + 8(p_g) 2 7 _(l T 7-)(| T 7- 
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Rearranging the various terms we can write 



/"(A) 



±4 



7+ + 7- 



(p + qYli(l T 7+) + {P-qYr-(l T 7-) 



(A.5) 



=F 



7+7- 



=F7+) =F7- 



2T7+T7- 



7+ + 7- 

(p+g) 2 7+(l T 7+) 2 + (P-g)V(l =F 7-) 2 
(1=F7+)(1=F7 



+ 



1 



=F 



+ 



2 =F 7+ =F 7- 



P+<?7+ + P-<?7- 



1 



2 =F 7+ =F 7- 7+ + 7 



2g 2 7+ (lT7+) + V7-(lT7-! 



The term in the last line is positive and bounded above by 4q 2 , both in 
the fermionic and bosonic case. For an upper bound in the fermionic case, 
we use that p±J± < D z , \p±\~/± < \fD~z as well as < j± < 1 to get 
/"(A) < 4g 2 (l + 2D Z ). Similarly we can obtain a lower bound. Using 
that p + q + p-q = 2pq in the second line in (|A.5|) . a simple estimate yields 
/"(A) > —12q 2 D z — 8p 2 q 2 in the fermionic case. Using these bounds in HA.lj) 
proves (|5.18j) . 

In the bosonic case, we only need to prove a lower bound on (|A.5|) . 
Proceeding as above, this time using p±7±(l + 7±) < D z and |p±|7± < y/Dl 
we obtain again /"(A) > —Ylq 2 D z — 8p 2 q 2 . This finishes the proof of the 
lemmas. ■ 
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